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Abstract 

We apply the method of QCD light-cone sum rules to calculate nonfactorizable contributions 
to the B — ► J/i/jK decay and estimate soft nonfactorizable corrections to the 02 parameter. The 
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oretical considerations and experimental data. Our result also confirms expectations that in the 
color-suppressed decay nonfactorizable corrections are sizable. 
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I. INTRODUCTION 



In nonleptonic decays of a B meson one can study effects of hadronization, perturbative 
as well as nonperturbative dynamics, final state interaction effects and CP violation. Mea- 
surements of exclusive nonleptonic B decays have reached sufficient precision to challenge 
our theoretical knowledge on such decays. It became clear that calculations have to reduce 
their theoretical uncertainties in order to make real use of data. Nowadays there exist sev- 
eral approaches which shed more light on the dynamical background of exclusive nonleptonic 
decays. The most exploited ones are Q CD factorization Q] and PQCD approach j^J. The 
PQCD model assumes that the two-body nonleptonic amplitude is perturbatively calculable 
if the Sudakov suppression is implemented to the calculation. In QCD factorization one 
can show the factorization of the weak decay amplitude at the leading order 1 / ra& level and 
can consider systematically perturbatively calculable nonleading terms of l/m& expansion. 
None of these approaches can take nonperturbative 0(l/m&) terms into account, but there 
is no evidence that such terms are negligible. 

The B — ► J/tpK decay is interesting because of the several reasons. There is a large 
discrepancy between the experiment and the (naive) factorization prediction. The naive 
factorization is based on the assumption that the nonleptonic amplitude (obtained in terms 
of matrix elements of four-quark operators by using the effective weak Hamiltonian) can 
be expressed as a product of matrix elements of two hadronic (bilinear) currents. It also 
predicts vanishing matrix elements of four-quark operators with the mismatch of the color 
indices. The naive factorization hypothesis has been confirmed experimentally only for 
class-I B -> D&M (M = it, p, a u D s , D*) decays. On the other hand, B -> J/tpK is 
the color-suppressed (class-II) decay and therefore a significant impact of nonfactorizable 
contributions is expected. 

Effects of a violation of the factorization hypothesis in the B —>■ J/ if)K mode have been, up 
to now, calculated by using different theoretical methods, resulting in the sign ambiguity of 
the decay amplitude i.e. 02 parameter (02 parameter is the effective coefficient of four-quark 
operators in the weak Hamiltonian; it is defined below by eqs. (|14|) and Q15|)). The QCD 
sum rule approach e 3j predicted a negative value for 02, while the PQCD hard scattering 
approach [^| and the calculation done in QCD factorization [5[ gave the positive value for 
the a2 parameter. Moreover, a detailed analysis of the experimentally determined B meson 
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branching ratios, although by assuming the universality of the a 2 parameter, gives conclusive 
evidence that generally the a 2 parameter should be positive 6]. On the contrary, the negative 
value of a 2 would indicate that the 1/N C term and the nonfactorizable part in the amplitude 
tend to cancel and would therefore confirm the large N c hypothesis Q. The validity of 
this hypothesis was established in two-body D meson decays, while, up to now, different 
attempts failed to prove this assumption for B decays (see i.e discussion in |8j). However, 
the sign ambiguity of a 2 cannot be solved experimentally by considering the B — > J/ipK 
decay alone. One of the possibilities is to consider the interference between the short- and 
long-distance contributions to B — > Kl + l~ 

Nonfactorizable corrections due to the exchange of hard gluons were calculated at 0(a s ) 
in QCD factorization. In this paper we concentrate on the LCSR estimation of soft nonfac- 
torizable contributions in the B — > J/ifiK decay coming from the exchange of soft gluons 
between the J/ip and the kaon. The calculation is based on the light-cone sum rules (LCSR) 
method [10]. This method enables a consistent calculation of nonperturbative corrections of 
hadronic amplitudes inside the same framework reducing therefore the model uncertainties. 

The paper is structured as follows. First, we discuss the results of the (naive) factorization 
method in Section |Hj The LCSR for the B — > J/ipK decay is derived in Section IIHI 
Next, to show the consistency of the method, we prove the factorization of the leading 
order contribution in Section IIV1 In Section the calculation of the soft nonfactorizable 
corrections is done by including twist-3 and twist-4 contributions. Section IVTI assembles the 
numerical results. In Section IVH| using the results of calculation, we discuss impact of the 
nonfactorizable term on the factorization assumption and the implications of the results. A 
conclusion is given in Section fVIIIt 

II. FACTORIZATION HYPOTHESIS AND NONFACTORIZABLE CONTRIBU- 
TIONS 

The part of the effective weak Hamiltonian relevant for the B — > J/ipK decay can be 
written in the form 

h w = ^v cb v; s [dMa + c 2 (fi)o 2 ] , (i) 
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which can be further expressed as 

G F 



where 



2 = (<a»(*r"&) ,6 2 = (cT^c)(sr^b) 



(2) 



(3) 



Here = 7 M (1 — 75), V^-'s are CKM matrix elements and A a 's are SU(3) color matrices. 
C*i(/i) and C 2 ([i) are short-distance Wilson coefficients computed at the renormalization 
scale fi ~ 0(m b ). The 2 operator appears after the projection of the color-mismatched 
quark fields in 0\ = (cT^b) (sT^c) to a color singlet state: 



O l = —0 2 + 20 2 . 



(4) 



The 1/N C term is the origin of the factor 1/3 in (J2J) . 

Under the assumption that the matrix element for the B — > J/i^K decay factorizes, the 
matrix element of the 2 operator vanishes because of the color conservation and the rest 
can written as 



(J/4>K\H W \B) = ^V cb V c : 
x(J/^K\0 2 \B) fact 



CM + 



CM 



l + 0(a s ) + 



A 



QCD 



(5) 



where the second and the third term represent hard and soft corrections to the factorizable 
amplitude, respectively. 

The factorized matrix element of the operator 2 is given by 



(J/^(p)K(q)\0 2 \B(p + q)y act = (J/^(p)\cT,c\0)(K(q)\m\B(p + q)) 



2e-qm J /^fj/ il F^ K {m 



J II" 



(6) 



where meson momenta are explicitly specified and p 2 = rn 2 ^. The J/ip decay constant is 
defined by the relation 

(J/^O)|c7 M c|0) = fj/^mj/^ (7) 

with e M being the J/ip polarization vector which satisfies the condition e ■ p — 0. The /ju 
denotes the J/ip decay constant determined by the experimental leptonic width T(J/ip — > 
l + l~) = 5.26 ± 0.37 keV by using the leading order calculation: 



/j/^ = 405±14MeV 



(8) 



The Fq K form factor is defined through the decomposition 



(K( q )\s llx b\B(p + q)) = (2q + p),F+ K (p 2 ) + ^ 2 P/j.(~^bk(p ) + ^bk(P )) (9) 



and estimated from the light-cone sum rules 
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12 1 has the value 



i^(m^)= 0.55 ±0.05. (10) 



By neglecting corrections in (J5|). the (naive) factorization expression for the B — > J/ipK 
decay emerges. Taking into account the NLO Wilson coefficients calculated in the NDR 
scheme Q for // = rn b {m b ) = 4.40 GeV and = 225 MeV, 

Ci(mF(m 6 )) = 1.082 C 2 {fr^{m b )) = -0.185 (11) 

and using the B meson lifetime r(5 ± ) = 1.653 ± 0.28 ps, we obtain for the branching ratio 
in the naive factorization 

B{B -> J/*fjK) fact = 3.3 • 10" 4 , (12) 

with the uncertainties in the order of 30%. This has to be compared with the recent mea- 
surements |3] 

B(B + -> J/i)K + ) = (10.1 ±0.3 ±0.5) • 10" 4 , 

B(B° -> J/ipK°) = (8.3 ± 0.4 ± 0.5) • 10" 4 . (13) 

Obviously there is a large discrepancy between the naive factorization prediction (JT^j) and 
the experiment. 

Returning to the expression for the B — ► J/ipK amplitude the corrections are given as 
an expansion in l/m& and a s . Apart from the 0(a s ) corrections to the factorizable part, there 
are also nonfactorizable corrections, which can be, either due to a hard gluon exchange or due 
to a soft gluon exchange (denoted in (j3J) as 0(a s ) or 0(Aq C D/ m b) corrections, respectively) 
between J/ijj and the B — K system. 

To be able to discuss the impact of the nonfactorizable terms, it is usual to parameterize 
the (J/i[)K\Hw\B) amplitude in terms of the a 2 parameter as 0,0] 

(J/if>K\H w \B) = V2G F V cb V; s e ■ qm m f m F+ K {m 2 j^) a 2 . (14) 
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The effective parameter a 2 is defined by 

a2 = g2M + giM + 2glM J^M . (15) 

The part proportional to the F^ K represents the nonfactorizable contribution from the 2 
operator 

(J/iPK\6 2 (n)\B) = 2e ■ qmj^f J/lp F+ K (fi 2 ) (16) 

and Fg K = corresponds to the naive factorization result, Eq. (jHJ). 

Because there is no explicit /x dependence of matrix element (jHJ), the /1 dependence of 02 
needs to be canceled by the nonfactorizable term. The nonvanishing nonfactorizable part 
is also required in order to suppress the strong renormalization scheme dependence of the 
effective parameter a 2 jl^ . 

Using the parameterization (|14p we can extract the a 2 coefficient from the experiment. 
From the measurements ()13j) one obtains 



exp I 



0.29 ±0.03, (17) 



with the undetermined sign of a 2 . 

On the other hand, with the NLO Wilson coefficients from the naive factorization 
yields 

4:nlo = C 2 {m b ) + <^p) = 0.176 , (18) 

which is significantly below the value extracted from the experiment. 

Following [gj], in Fig.l we show the partial width for B — > J/ipK as a function of the 
nonfactorizable amplitude Fbk- The zero value of Fbk corresponds to the factorizable 
prediction. There exist two ways to satisfy the experimental demands on the Fbk- According 
to the large 1/N C rule assumption [7], one can argue that there is a cancellation between 
1/N C piece of the factorizable part and the nonfactorizable contribution, ()15|) . That would 
demand relatively small and negative value of Fbk- The other possibility is to have even 
smaller, but positive values for Fbk, which then compensate the overall smallness of the 
factorizable part and bring the theoretical estimation of a 2 in agreement with experiment. 

One can note significant /i dependence of the theoretical expectation for the partial width 
in Fig.l, which brings an uncertainty in the prediction for Fbk{h) of the order of 30%. This 



6 



FIG. 1: The partial width T(B — > J/ipK) as a function of the nonfactorizable amplitude Fbk- 

uncertainty is even more pronounced for the positive solutions of i<W(/x). The values for 
Fg K extracted from experiments 

F+ K (m b ) = 0.028 or F+ K (m b ) = -0.120 , 
F£ K (m b /2) = 0.046 or F+^(m 6 /2) = -0.095 . (19) 

clearly illustrate the \x sensitivity of the nonfactorizable part. 

In the following we will calculate the nonfactorizable contribution F^ K , which appears 
due to the exchange of soft gluons, by using the QCD light-cone sum rule method. 



III. LIGHT-CONE SUM RULE FOR (J/ipK\0\B) 



A. The correlator 



To estimate the soft-gluon exchange contributions to B — > J/ibK we use the method 

n 

developed in [10J for the B — > txtx case. In this approach one considers the correlation 
function: 

F v (p, q, k) = i 2 [ d'xe-^> [ C /V (p -^(^(g)|T{j^(y)O(0)jf (x)}\0) , (20) 



where j J J^ = c^ v c and jf = m b bi^j^u are currents interpolating the J ftp and B~ meson 
fields, respectively. The correlator is a function of three independent momenta, chosen by 
convenience to be q, p — k and k. Diagrammatically the correlator is shown in Fig.2. 
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K(q) 




FIG. 2: B — > J/ipK decay in LCSR. The shaded oval region denotes nonperturbative input, the 
K meson light-cone distribution amplitude. J ftp and the B meson are represented by the currents 
j J /^(j) — k) and j B {p + q), respectively. The square stands for the Oi four-quark weak operators. 

Here, it is important to emphasize the role of the unphysical k momentum in the weak 
vertex. It was introduced in order to avoid that the B meson four-momenta before (ps) 
and after the decay (P) are the same, Fig.2. In such a way, one avoids a continuum of light 
contributions in the dispersion relation in the 5-channel. These contributions, like DD* S 
or D*D S , have masses much smaller than the ground state B meson mass and spoil the 
extraction of the physical B state. Also, they are not exponentially suppressed by the Borel 
transformation (see for example the discussion in 

The correlator (j2*Uj) for nonvanishing k is a function of 6 independent kinematical invari- 
ants. Four of them are taken to be the external momenta squared: (p + q) 2 , (p — k) 2 , q 2 
and k 2 , and additionally we take P 2 = (p — k + q) 2 and p 2 . We neglect the small corrections 
of the order 0(m 2 K /m 2 B ) and take q 2 = m? K = 0. Also, we set k 2 = 0. The p 2 momentum 
is for the moment kept undefined, in order to be able to make unrestricted derivation of 
the sum rules. Its value is going to be set later, by considering the twist-2 calculation of 
the factorizable part, Section ITVl and will be chosen p 2 = fn 2 j/^ in order to reproduce the 
factorization result, (jUJ). Furthermore, we take (p — k) 2 , (p + q) 2 and P 2 spacelike and large 
in order to stay away from the hadronic thresholds in both, the J / if) and the B channel. All 
together we have 

q 2 = k 2 = 0, p 2 = m 2 !^ \{p-k)\ 2 ^> A Q cd, \(p + q)\ 2 > Aqcd, \P\ 2 > Aqcd • (21) 
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B. Derivation 



The first step is the derivation of the dispersion relation from the correlator (|2U|). Inserting 
a complete set of hadronic states with the J / ip quantum numbers between the J/ tp current 
and the weak operator in (|2Up gives us the following: 



E< (A)n ((p+^) 2 ' p2 ^ 2 ' e(A) -* 



J a HJ/w s-(p-k) 2 

where qi — p,q,k and e^ x ' ■ (p — k) = 0. The sum runs over the polarizations of J/ip. The 
lowest state contribution satisfies 

(0|757„c|J/V(p - k, e^)) = m m fj,^ (23) 

and (p — k) 2 = m 2 ,^. In ()22jl. and Sq'^' are the spectral density and the threshold 
mass squared of the lowest excited resonances and continuum states of the J/ip channel, 
respectively. 

The hadronic matrix element of interest is denoted by 

U((p + q) 2 ,P 2 ,p 2 ,e^-q l )= i J d A xe-^ + ^(J/^{p-Ke^)K{q)\T{0^)ji{x)m. (24) 

On the other hand, for spacelike (p — k) 2 <C m j/ 1 / ) ^ ar away from the poles associated with 
the resonances and continuum states, the correlator F v can be calculated in QCD in terms 
of the quark and gluon degrees of freedom and written in a form of a dispersion relation as: 

Tu = \T i; > s .^( g , (p + q) 2 ,P 2 ,p 2 ) 
7i J 4m 2 s-(p-k) 2 

with the kinematical decomposition 

T v = (p - k) v F^ + k u F^ + q u F® + e vaPl {p - k) a kfq^F^ . (26) 

By assuming quark-hadron duality one substitutes the hadronic spectral density p\[^ in 
(J22J) with the one calculable in QCD and replaces s^ J ^^ with the effective threshold of the 
perturbative continuum, Sq , i.e : 

ptfis, (P + tf, P\ P 2 M* - = -l^sMs, (P + ?) 2 , P 2 , P 2 Ms - 4^) . (27) 

7T 

9 



By matching the hadronic relation (|22|) with the QCD calculation (|25|) one obtains the 
sum rule expression 



mj/4>fj/ip 
m JM - (P - fc ) 2 



1 ^ Jm s ^(s,(p + g) 2 ,P 2 ,p 2 ) 

7T 7 4m 2 S - (p - k) 2 

In order to reduce the impact of the approximation ()27|) and to suppress contributions 
from excited and continuum states, as usually done for quarkonium systems one performs 
n derivations in the momentum (p — k) 2 and receives n- moment sum rule for the correlator 
n((p + q) 2 , P 2 ,P 2 , ■ Qi) of the form 



£< {A) n((p+^pv,e (A) ^. 



J J^ (™2 i r,2\n+l 



■V-T rfg ( 7 7 t nSli^ (P + P 2 ) . ( 29 ) 

where Qo is the sum rule parameter that role will be discussed later, in Section IVT1 

We proceed by using the analytical properties of the II((p + q) 2 , P 2 ,p 2 , • qi) amplitude 
in the (p + q) 2 variable of the P-channel and insert in (|24|) the complete set of hadronic 
states with the B meson quantum numbers which yields 

E< (A)n ((p+^ p2 ^ 2 ' e(A) -*) = 

A 

2 V2 J2e^(m P -Ke^)K(q)\O(0)\B(p + q)}+ r ^^'f, 1 . 

m%-{p + q) 2 ^ J s kb) s' - (p + g) 2 

(30) 

In above, as before, it is assumed that in the last term the polarization sum is already done. 

The QCD part, given by the r.h.s of the eq. (|29)l and rewritten in a form of the dispersion 
relation, now in the (p + q) 2 variable, exposes the form of the double dispersion relation as 

1 ft" . {rrfju + Ql) n+1 



r j ° f \ (m 2 J/ ,+Q 2 T + \ T( ( ^ , 2 p2 2 , 

/ ds (TTT^l l ™sfu(s, (p + q) 2 , P 2 , p 2 ) 

i f 4 '* , K« + Q 2 X + ' r'^ p2 ^ ds' 



From the Maldestam representation of the kinetic variables one can see that the integration 
limit of s' variable is going in general to depend on s, P 2 and p 2 and we denoted these depen- 
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dence by /i(s, P 2 ,p 2 ). In the following, those terms which disappear after taking moments 
in J I ip channel and after making the Borel transform in the i?-channel are neglected. 

In order to subtract the continuum of B states, we exchange the order of the integration 
in (|31|) and use quark-hadron duality in B channel in a sense that the spectral density pf 
is approximated by the s' > s B part of the double dispersion integral (j3Tj) . where s B is the 
effective threshold of the perturbative continuum in the B channel. Therefore, 



m 



ffl r 2 E < {X) ( J /^(P - ^ e^)K(q)\O(0)\B( P + q)) 



1 



J/i> i 2 



+ Ql) n+1 rM-tfS*J) ds' 



x/ da \ J/ *ZL / - ; 8 , Tm.ImXi^^ 2 ,/) (32) 



Uml (s + Qo)" +1 Jml S'-( P +q) 2 

and after the Borel transformation in (p + q) 2 variable, we can further write 

1 



E < (X) ( J /^(P - k > e^)K(q)\O(0)\B(p + q)) 

A 



x P f rfs 'e«-.'>/-I m „I m ,^( S , P^) , (33) 

J Ami \ s + VoJ ./mg 

In above, M is the Borel parameter and the function f 2 is the upper limit of the s integral 
after subtraction of continuum of B channel. 

Further, to extract the kinematical structure of interests, we decompose the matrix ele- 
ment {J/ip{p-k,eW)K(q)\O(0)\B(p + q)} as 

(J/i/)(p- k,e (x) )K{q)\O(0)\B(p + q)) = e • q A® + e • k A (k) 

+e a e ap ^{p-k) p k a q^A ie) (34) 

By inserting this expansion in the expression (}3*3*|) . after the summation of the polarization 

vectors 

C *(A) J\) _ ( _ n , (p-k)v{p-k) c 



Z&w = (-^+ toZr > (35) 

A x 

one obtains the sum rule for different kinematical structures: 



-k v A™ - q v A® - e^(p - k)"(f0 A® 

+(p - *o* ( ip ~ k) ,;, k ^ + ip r k) : q AW 



\ (p — A;) 2 (p — A;) 2 

1 r'°* (m 2 T/h + Ql) n+1 rMs,s^PV) /2 
/ ds—^————— / ds'e {m B~ s)IM 

xIm s Jm s [k u + q v + (p - k) v F^ + e vpai {p - k) p q a k^ F®] 

(36) 
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The coefficient function in front of (p — k) v looks like 

A ( P -k) _ ( (p-k)-k (k) (p-k)-q ( ) \ 

A -\(P-ky A + (p-ky ) ' (37) 

which is a consequence of the conserved J/ip current. The sum rule expression for the A^ - ^ 
part reads 

1 



j\(.p~ k ) 



B 



dl* (Vn 2 i Q2\n+l f2{ss Bp2 2) 

x / ds { f ^K-^ Ims , Im8F M) (S)S ' )P 2 )p2)i 



4m? 



(s + Ql) n+l 



(38) 



At the end we analytically continue P 2 to P 2 > 0, and choose P 2 = m 2 B . This enables the 
extraction of the physical matrix element, because the unphysical momentum k disappears 
from the ground state contribution, due to the simultaneous conditions applied, P 2 = m 2 B 
and (p + q) 2 = m 2 B . From and (|37|) follows 

(J/iP(p,e^)K(q)\O(0)\B(p+q)) = e-qA^(P 2 = m|) = e-q ^ - A^ k \P 2 = m B ) (39) 

171 B ~ P 

and the final sum rule relation for the physical matrix element (J/ip(p, e^)K(q)\O(0)\B(p + 
q)) takes the form 

(J/1>(p,eW)K(q)\0(p)\B(p + q)) = 



(40) 



X 



^m a ,'hs\ 8 F {p - k) {s,s\m 2 B ,p 2 



2/2 2 
m J/i,{ m B ~P 



Some comments are in order. The B — > J / r xpK case seems to be much more complicated 
than the decay of a i? meson to two light pions discussed in ^(|. The complication does 
not appear only due to the massive c quarks, or the vector structure of the J ftp current, but 
mainly due to the local duality assumption in J / ip channel, which is expected to work much 
worse than in the pion channel in the B — > tttt decay. Although it is possible to stay away 
from the the excited and resonant hadronic states in the J/ip channel, one can still expect 
that there will be an influence of the ip' resonance, which, in a more precise calculation has 
to be taken into account explicitly. The technical difficulties which are induced by the fact 
that the value of P 2 parameter is close to the hadronic threshold of J/ip are left for the 
discussion in Section IVl 
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IV. FACTORIZATION IN THE LIGHT-CONE SUM RULE APPROACH 



We first consider the contribution of the O2 operator. As we have shown in the 
introduction, this operator contributes to the factorizable part of the matrix element 
(J/i}(p,eW)K(q)\H w \B(p + q)). 

The main contribution comes from the diagram shown in Fig.2, where for O = O2 there is 
no interaction between the charm loop and the B — K system at the leading level. Therefore, 
the calculation of this contribution is rather simple. According to the expression ([40)1 . the 
(p — k) v part of the correlation function ()20|) . F( p ~ k ^ needs to be calculated and its double 
imaginary part has to be extracted. The calculation proceeds in several steps. One inserts 
first explicitly the J/ip and B currents in (|2Uj) . and takes the expression (j3J) for the operator 
02- The c-quarks are contracted to a cc- loop and can be then independently integrated. 
The contraction of 6-fields produces a free 6-quark propagator and the rest of the fields is 
organized into the leading, twist-2 kaon distribution amplitude <pK- Explicitly, we obtain 



,(p-fc) _m 2 J K f°° ds ^ ^ ^ ^ t 2m 2 c \ ^ 4jn% f 1 ^ <j) K (u) 



FZ2=-r^ 7 j^q-ip-k) 1 + -^ W1 / du ., ( 

4vr 2 J 4m 2S-(p-k) 2 ^ J \ s J \ s J m 2 b -(p + quy 

(41) 

where 4>k(u) is the kaon twist-2 distribution amplitude defined by 

(K{q)\s{0) W u{x)\0) = -i % f K f due luqx Mu) ■ (42) 

Jo 

The first integral in ([41 p. apart from the kinematical factor q ■ (p — k), is nothing else but 
he charm loop contribution to the vacuum polarization calculated in the sum rule approach 
The second integral, considered in the leading twist approximation, reduces exactly 
to the light-cone twist-2 expression for the F^ K form factor Q|. This part, with the 
substitution u = (m 2 — p 2 ) / (s' — p 2 ), can be rewritten in a dispersion form as 

f 1 , <M«) r ds' 1 ± . . A , 

X dn ,nl-( P + q uY = J ml ,- {v + q y,-pM^ » ■ («) 

In such a way the expression (|41j) receives the needed double dispersion form from which 
the double imaginary part in s and s' variables can be trivially extracted. 

The contribution of the O2 operator to the B — > JipK matrix element then follows from 
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the the sum rule relation ()4Uj) : 



= 2e-q 



B{p + q)) 



tw2 



(m 



JH> 



+ Q 



2\n+l 

o) 



Air 2 



ds 



X 



fKm 2 h 



4m'} 



{* + Ql) n+l 



1 + 



2mt 



Am 2 p 



m 



B 



s m 



m (i _ p.) 

m B 



ds' 



s')/M 2 



2e- qmj/ffj/^F^x 



2 S' — p 2 



1 1 

s' — p 2 



(P 2 ) 



(44) 



Here we see that the amplitude (J/ip(p, e^)K(q)\O2(0)\B(p + q)) factorizes and in a good 
approximation the factorizable expression for the 

(J/ip(p,e^)K(q)\O 2 {0)\B(p + q)) amplitude, given by (0) , is recovered for p 2 
the first parenthesis, apart from the small (1 — s/m 2 B )/(l — m 2 j^/m 
the leading order expression for the f 2 ^ in the QCD sum rule approach. The correction 
is the result of calculation with the nonvanishing k momentum. The second parenthesis in 
()44|) gives the twist-2 contribution to F BK (p 2 ) form factor. By reproducing the factorization 



2 B ) correction, there is 



result for p 2 



m 



, we fix the value p 2 



l2 j/i> 



m%,.,. also in the further calculation. 



V. SOFT NONFACTORIZABLE CONTRIBUTIONS IN THE LCSR APPROACH 

For a discussion of nonfactorizable contributions to the B — > J/ ipK decay, we need to do 
a systematic a s and twist expansion of the correlator ()2())) . 

After explicit insertion of the interpolating J/ ip and B meson currents and the operator 
O2 or O2, the correlation function (J2*U|) can be written in form: 

(K{q)\Tr[ lu S % i(y, 0|m c )r M S fe4 (0, y\m c )\ s l T»S mn (0, x\m b ) l5 u n \0) , (45) 

where i, j, k, I, m, n are color indices, S(x,y\m) are quark propagators defined in (j4fij) below 
and = 5^ or = T* J = (X a /2)^ for the insertion of O2 or O2 operator, respectively. 

The a s and twist expansion is achieved by considering the light-cone expression for quark 
propagators. Up to terms proportional to G, the propagation of a massive quark in the 
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K(q) 



K(q) 




Jv ¥ (p-k) 



js(p+q) 




jf(P-k) 



FIG. 3: Soft nonfactorizable contributions to the correlation function 1)20(1 . 
external gluon field in the Fock-Schwinger gauge is given by 

S ij {x 1 ,x 2 \m) = -t(0\T{q i (x 1 )q j (x 2 )}\0) 
d 4 k 



(2tt) 



-ik(x\—X2) _ 



x 



ft + m 



2 (/c 2 — m 2 ) 2 Ml/ fc 2 — m 2 



V(X 1 - X 2 )^u 



(46) 



From the above, considering the color structure of the 2 operator, we can easily deduce 
that the nonfactorizable contribution from this operator appears first at the two-gluon level 
and is therefore of 0(a 2 ). On the contrary, nonfactorizable corrections from the (D 2 operator 
are already given by the one-gluon exchange. The leading hard nonfactorizable contributions 
are due to the exchange of a hard gluon between the c-quark (antiquark) and the one of the 
remaining b, u or s quarks, see Fig.2. These contributions emerge at the two-loop level and 
although they are calculable in LCSR, their calculation is technically very demanding and 
will not be discussed in this paper. 

Insertion of the gluonic term of the propagator (y,0,m c ) or S kt (0,y,m c ) yields the 
contributions represented in Fig. 3. These are the leading soft nonfactorizable contributions. 
In terms of the light-cone expansion they are of the higher twist and described by the three 
particle kaon distribution amplitudes defined by the following matrix elements: 
- twist- 3 distribution amplitude 



(Q\s(0)a^ 5 G af 3(vy)u(x)\K + (q)) = if 3 K [(qaq^gpu - q^g a u) 



(47) 
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twist-4 distribution amplitudes 

(0 MO)/-/;,, <('7/)»!-'-) K ' (>/)) ^ Ik I Va4^^)e-^' + ^ " 

qx 

-iq(xai+yvas) 



+(g m Qp- 9tfQa) / Vai<j)±(ai, fi)e 



(48) 



(0\s(0)l^G a p(vy)u(x)\K + (q)) = g M 



q a xp - qpx a 
qx 

iq(xai+yvag) 



f K / Va^a^^e-^^^ 



(49) 



zk{oci, fx) = 360aia2«3 ^1 + a(//)-(7a 3 - 3) + b(fj)(2 - 4ai« 2 - 8a 3 (l - a 3 )) 
+ c(/i)(3«ia 2 - 2a 3 + 3a 2 )) , 



In above, G a p = ^e a p pa G pa , G pa = g s X a /2G pa , Va>i = daida 2 da 3 S(l — «i — ct 2 — a 3 ), and 
9aB = 9aB ~ (x a q/3 + xpq a )/qx. Both twist-3 and twist-4 distribution amplitudes contribute 
at the same order. They are parameterized by 



(50) 
(51) 
(52) 
(53) 

(54) 

s. 



011 (aj,/x) = 1205 2 (/i)e(/i)(«i - a 2 )aia 2 a3 



- + 2e(/i)(l-2a 3 ) 



4>±((*i,fj) = 305 2 (/i)a 2 (l - a 3 ) 
0||(tti,/i) = -1205 2 (/i)aiQ; 2 Q! 3 



+ 2e(/i)(l-2a 3 
+ e(/i)(l - 3« 3 
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2 it ] and the values are listed in 



The parameters are estimated from sum rules 
In the numerical evaluation we use the asymptotic form of the above expressions where 
a(/i), b(fjL), c(/i) and e(/x) dependence is neglected. The asymptotic expressions for twist-3 
and twist-4 distribution amplitudes should provide sufficiently reliable estimates of already 
subleading contributions. 

The QCD calculation of two diagrams in Fig. 3 at the twist 3 level yields 



7 {p-k) 
tw3 



m b f: 



3K 



4vr 2 



dv I Vcti 



dx- 



2x 2 (l - x) 



m 2 b -{p + q(l - ai )) 2 J m 2 _ Q 2 x (i _ x ) 



x q ■ (p — k) (2 — v)q ■ k + 2(1 — v)q ■ (p — k) 



(55) 



where Q = p — k + va 3 q. Comparing the above expression with the one obtained for the 
B — > 7T7T case, (Eq. (26) in Q|), we can see that there is an additional, x-integral for the 
massive cc loop. Otherwise, the expressions are the same and for m c — > the result form 
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Eq. (26) in [.llj is exactly recovered, up to a sign, which can be traced back to a difference 
between the pseudoscalar and vector currents interpolating tt and J/ip, respectively. 

By changing the order and variables of integration one can bring the above expression 
into the following form: 



tw3 



m b f: 



3A 



16vr 2 



ds 



Ami 



dy 



s - (p - k) 2 J 2^/y y :c ( Sj2/ ,p2 ) m 2 b -(p + qu) 



du 



x 



^0 3 a(1 ~U,U-V,V) 

x(s,y,P 2 ) V 

» - — + ((P + ?) 2 " P 2 ) (2« - x(s, y, P 2 )) 



x 



and 



1-2/ 



x(s,y,P 2 



(56) 



4m 2 

T=7 



.-!» • (57) 
It is important to emphasize here that the above expression (|56j) is defined only for large 
spacelike momentum |P 2 | ~ m 2 . Furthermore, the expression (JSfiJl does not have a needed 
double dispersion form. 

In order to proceed we write 



f/.S' 



2 d2 



^ itf S ~ (p - /C) 2 



(58) 



where 



s,(p + g) 2 ,P 2 ,p 2 ^ 



3A 



167T 



X 



(i 2^ J x (s,y,p2) m 2 b -(p + qu) 2 J x ( s ,y,PZ) v 

Am 2 



^-03a(1 -u,u-v,v) 



s — 



y 



+ {(p + q) 2 -p 2 ) (2v-x(s, yj P 2 )) 



(59) 



Now, it is possible to use the quark-hadron duality in J/if> channel and to subtract the 
J /if) continuum states by approximating them by (}59|) . which changes the upper limit of 



s integration in (JH%|) to s J q^ ~ 15GeV 2 . This restriction of the s integration enables the 
expansion of the imaginary part lm. s F tW 3 in the x(s, y, P 2 ) variable. To reach the satisfactory 
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precision we expand (JHHJ) up to order 0[ 



x 



lm s F t i; k \s,(p + q) 2 ,P 2 



du 




o m? b - (p + qu) 2 J 
3K (1 -u,u-v,v) 
u dv 



- _dy_ 

Ami 



+ 2v((p + q) 2 -p 2 



—(j)3 K {l -u,u-v,v) (0 + q) 2 - p 2 ) 



Am 2 



b 3K (l -u,u-v,v t 



v=0. 



Am 2 



d 



x(s,y,P 2 ) 
x 2 (s,y,P 2 ^ 



+ 0(x 3 



(60) 



v=0 



In the above expression it is important to keep in mind that s receives values in the 
range Ami < s < Sq . It has also to be noted that the coefficients in the expansion are 
P 2 independent objects. So, although, the above expression was derived for P 2 < 0, the 
complete expression (|4U|) for the physical amplitude B — > J/ipK is an analytic function in 
P 2 and it can be analytically continued to the positive values of P 2 = m 2 B . The result 
is more reliable for smaller 0(s/P 2 ) corrections. In our case, although the expansion is 
well converging, the first order correction in ~ x(s,y,P 2 ) amounts to ~ 25%, which is 
significantly larger than the similar correction ~ Sq/P 2 ~ 0(lGeV/m 2 B ) in the B — > 



case. Therefore, in the calculation of the soft nonfactorizable correction for B — > J/ipK, 
the analytical continuation of P 2 to its positive value embeds an unavoidable theoretical 
uncertainty. However, 0(x 2 ) corrections are already at a percent level, and the expansion is 
well converging. 

The same procedure employed for twist-4 contributions gives somewhat more complicated 
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result: 



,(p-k) _ ™lf. 



K 



ds 



dy 



J iml s - (p - k) 2 J 2^/y J x{S) y ]P ^ m 2 b -(p + qu) 



du 



mlf K 



— U,U — V,V) 

(s,y,P 2 ) V 



3 x(s,y,P 

v 



ds 



lir 2 J 4m 2 s-(p-k) 2 

-5-^1(1 -u,v) 

1* 



x(s,y,P 2 ) v 



m 2 b f K 



8tt 2 



f/.S 



4m^ 



dy 

2 Vy Jx(s,y,P 2 ) [(ml - (p + g«) 2 ' 
s - y-^ + ((p + g) 2 - P 2 )(-v + x(s, y, P 2 )) 

1 dy * 



du 



x- 



it- 



s - (p - A;) 2 7 2^ J x { s ,y,p*) [m 2 - (p + gw) 2 ] 

+ ((P + # " P 2 ) (-« + x(s, y, P 2 ) 



i-y 

ds 



m 2 b f K 

' 8vr 2 i 4m§ ^-(p-A;)^ 



dy l-y 
2 Vvi-y-^ 



du 



x(s,y,P*) m l ~ (p-quf 



u 



x(s,y,P 2 
-P 2 



(2q-{p-k)) 2 (l- 



x(s, y, P 2 ) q ■ k 
u q ■ p 



(61) 



Here, q ■ k = |((p - k) 2 - P 2 + (p + q) 2 - p 2 ) and q ■ p = |((p + q) 2 - p 2 )). 
The twist-4 wave functions appear in combinations 



$i(tt,i>) = / do; ( 4>±(ui, 1 — a; — t> , t> ) + <j>\\(u), 1 — — v, v) 
$ 2 (w) = 



da/ 



1-uj' 



du" </>±(w", 1 - u" - u', 00') + 0|| (w", l-u" - u', 00') (62) 



The first term in (|6T|) can be treated in a similar way as the twist-3 part, F^ 3 k \ expanding 
in x(s,y,P 2 ) with the result 



i(p-*) 



A' 



tu)4 



ds 



7T 2 74m? s - (p - A;) 2 7 m 2 - (p + qu) 



du 



dy 

2 Vv 



l Udv l n ^ 
x / — 0±(1 - u, u - v,v) 



3 x{s,y,P 2 



+ 0(x 3 ). 



(63) 



Other parts in (JfiT|) contain denominators of a form 

1 1 



[*-(p-A;) 2 ] 2 
which are typical for twist-4 contributions 



or 



[ml - (p + uq)'- 



i2 ' 



(64) 
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To be able to deal with such terms we perform a partial integration. However, the problem 
is the subtraction of a continuum for such terms, because the complete expression does not 
possess the needed dispersion form, where the hadronic spectral density can be identified 
with the imaginary part of the QCD amplitude, unless the surface terms are equal to zero. 
Fortunately, twist-4 contributions with the higher power of denominators numerically appear 
to be suppressed. Their contribution, neglecting the surface terms, is in the region of a few 
percent. Uncertainties involved in the LCSR calculation are certainly much larger, and we 
argue that the contributions with the higher power of denominators in the twist-4 part can 
be safely neglected in the numerical calculation. 

It is important to emphasize that due to the specific configuration of momenta, imposed 
by the J/ip continuum subtraction, the analytical continuation of P 2 does not produce an 
imaginary phase. From this continuation, one would expect to get an imaginary phase in 
the penguin contributions of operators 0\ and C 2 . The phase is typical for such kind of 
contributions and known as the BSS phase Q]. However, the penguin contributions in the 
process under the consideration are suppressed in the large N c limit by 1/N C (additionally 
to the l/N c suppression of the emission amplitude calculated here), and are beyond a scope 
of this calculation. 

Putting twist-3, eqs.()58|) and (pOjl. and twist-4 (J53)l expressions together and subtracting 
the continuum of B states, the final expression for the soft contributions to the B —>■ J / ipK 
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amplitude, in the approximations discussed above, has the form 



(J/i>(p,eM)K{ q )\O{0)\B{p + q)) = 2e-qm m fj M 



1 



[m j/^ + QlT +1 1 f 1 du (m 2 B -(m 2 b -m 2 J/4> (l-u))/u)/M 2 



ds- 



^ 2 f 2 j,J± m 2 (s + Q 2 o) n+1 2m%f B J u B 
« dy m b 

%Vy m B - m2 ,/ip 



u 



3A' 



"* 1 
— 



v- 



m b ~ m2 J/ib 

1 — u, u — v, v) I — (2t> — x(s, y, m B )) + s 



ii 



- 18- 



Am 2 c 



1 



2 r- 3 i<r(l - u, u - v , v) ) x(s,y,m 2 B ) 
1 — y J \tr 



+m b f J 



Ami 
u dv~ 



dv \ v 2 



v=0 



> 3K (1 -u,u-v,v) 



x 2 (s,y,m 2 B ) 



J v=Q 



K 



'1 — u.u — v, v) 



3 - -x(s,y,m 2 B ) 



where u% = (m 2 - m^)/(s£ - m 2 M ). 



Am c 



(65) 



VI. NUMERICAL PREDICTIONS 



Before giving numerical predictions on the soft nonfactorizable contributions, we have 
first to specify the numerical values of the parameters used. 

For parameters in the B channel we use mg = 5.1 GeV and the values taken from jll|: 
f B = 180 ± 30 GeV, m b = 4.7 ± 0.1 GeV, and s% = 35 ± 2 GeV 2 . For J/ if) we use the 
following: mj^ = 3.5 GeV, fj/^ = 0.405 ± 0.014 GeV from Eq.flEJ), m c = 1.25 ± 0.1 and 
s J q^ = 15 ± 2 GeV 2 [l^. The K meson decay constant is taken as fx = 0.16 GeV. For 
parameters which enter the coefficients of the twist-3 and twist-4 kaon wave functions we 
suppose that / 3vr ~ f 3K and 5 2 K ~ <5 2 , and take f 3K = 0.0026 GeV, 5 2 (fi b ) = 0.17 GeV, where 



m 2 ~ 



m b /2 ~ 2.4 GeV |2(j,|21|. 



Like in any sum rule calculation it is important that the stability criteria for (|4"UJl are 
established by finding the window in n and M 2 parameters in which, on the one hand, excited 
and continuum states are suppressed and on the other hand, a reliable perturbative QCD 
calculation is possible. The stability region for the Borel parameter is found in the interval 
M 2 = 10 ± 2 GeV 2 , known also from the other LCSR calculation of B meson properties. 
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Concerning moments in J/ip channel, the calculation is rather stable on the change of n in 
the interval n = 4 — 6. Qq is parameterized by Qq = 4m^, where £ is usually allowed to take 



values from to 1. As it was argued in 17(, where sum rules was applied for calculating the 
mass of J/ip, and was also observed in our calculation, at Qq = (£ = 0) there is essentially 
no stability plateau where n is small enough that the QCD result is reliable and at the 
same time the lowest lying resonance dominates. More stable result is achieved for £ ^ 0. 



However, the result appears to be sensitive at most to the variation of the parameters Sq 
and s J q^ . 

The numerical results for the soft nonfactorizable contributions are as follows 

*BK,*M = °- 0051 > *$K,*M = °- 0089 > (66) 

calculated at the typical /x& ~ rrib/2 scale of LCSR calculation. The above values are 
obtained for n — 5, M\ = 10 GeV 2 and £ = 0.5. In general, one could expect that twist- 
4 contributions are relatively 0(l/m&) suppressed with respect to the twist-3 part and 
therefore are smaller. However, careful study of the heavy-quark mass behavior of the final 
expression (|65|) shows that in the heavy-quark limit the twist-3 and twist-4 contributions 
are of the same order (and are both suppressed by l/m& with respect to the factorizable 
part (|44jl ) . Therefore, it is not surprising that the numerical contribution of the twist- 
4 part is relatively large. Even- and odd-twist contributions stem from different chiral 
structures of the 6-quark propagator and are, therefore, independent. The l/m& suppression 
should, however, certainly be true when we compare even(odd)-twist contributions among 
themselves (i.e. the twist-4 with the twist-2 contribution; the twist-5 with the twist-3 part 
etc. ). 

The variation of the sum rule parameters implies the values: 

*BK,tu*W = °- 004 - °- 007 > FBK,tM = 0-006 - 0.012 (67) 
and the final value 

F+ K (fx b ) = 0.011 - 0.018 . (68) 

First, we note that the nonfactorizable part (|68|) is much smaller than the B—K transition 
form factor (fTU|) which enters the factorization result. It is also significantly smaller than its 
value (|TT?|) extracted from experiments. Nevertheless, its influence on the final prediction for 
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a2 is significant, because of the large coefficient 2C\ multiplying it. Furthermore, one has to 
emphasize that Fg K is a positive quantity. Therefore, we do not find a theoretical support for 
the large N c limit assumption discussed in Section HI that the factorizable part proportional 
to C\(fi)/3 should at least be partially canceled by the nonfactorizable part. Our result 
also contradicts the result of the earlier application of QCD sum rules to B — > J/ipK jjjj], 
where negative and somewhat larger value for Fg K was found. However, earlier applications 
of QCD sum rules to exclusive B decays exhibit some deficiencies discussed in In jsj, 
mainly the problem was the separation of the ground state contribution in the B-channel 
and the wrong m& — > oo limit of higher-twist terms obtained by using the short- distance 
expansion of the four-point correlation function. In this work, following the procedure taken 
from ^(|, the problem is solved by introducing the auxiliarly momentum k in the 6-decay 
vertex and by applying the QCD light-cone sum rules. 

Using the same values for the NLO Wilson coefficients as in Section |Hj one gets from 
(|68|l for the effective coefficient a 2 the following value 

a 2 ~ 0.15 - 0.18 . (69) 

Although the soft correction contribute at the order of ~ 30% — 70%, the net result (fBTJj) 
is still by approximately factor of two smaller than the experimentally determined value 



(J/i>K\(D\B) = (K\sTMB)(J/il;\cr v c\0) 



1 



VII. QCD FACTORIZATION FOR THE B -» J/ipK DECAYS AND THE IMPACT 
OF SOFT NONFACTORIZABLE CORRECTIONS 

In an expansion in l/ra& and a s , matrix elements for some of two-bodv decays of a B 
meson can be computed consistently by the QCD factorization method [l|. This model 
applied to the B — > J/ipK decay gives 

+ 0(a s ) + o(^ 

V m b 

= F BK (m 2 j^) I T 1 (u)(j)j/^(u) 
Jo 

+ J didudvT 11 ^ u, v)<f> B (0<f>K(v)<f>J/4 u ) + (^p) • ( 7 °) 

T 1 and T 11 are perturbatively calculable hard scattering kernels and 4>m=b,k,j/tt are meson 
light-cone distribution amplitudes. T 1 starts at order 0(a°), and at higher order of a s 
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contains nonfactorizable corrections from hard gluon exchange or penguin topologies. Hard 
nonfactorizable corrections in which the spectator of B meson contribute are isolated in T 11 . 
Soft nonfactorizable corrections denoted above as 0(AQCD/ m b) effects cannot be calculated 
in the QCD factorization approach. According to some general considerations Q these 
effects are expected to be suppressed, but there is no real confirmation of this conclusion. 

In the limit m c ~ mj — > oo, it can be shown that at the leading order in l/my, there is 
no long distance interactions between J / ip and the rest B — K system and the factorization 
holds. Actually, the J/ip case is somewhat exceptional, since soft gluons in this limit are 
suppressed only by a factor Aqcd/(^6«s) 111 rather than by Agc^/m?, like, for example, in 
the B — > Dtt decay, for which the factorization has be proved at the two- loop level. If J/ip 
is treated as a light meson relative to B, then the factorization is recovered at m c /m/, — > 
limit. Unfortunately, for the higher l/m?, corrections, the factorization breaks down 

In connection to (J7UJ) the following should be emphasized. In the heavy quark limit 
mj — > oo the hard scattering kernel T 1 is nothing else but the J / if) meson decay constant and 
by neglecting a s and 0(AQCD/mb) corrections, the naive factorization result (0) is recovered. 
In the hard corrections appear J/ if) and B meson light-cone distribution amplitudes. Under 
the assumption that m c <C mj, the light-cone distribution amplitudes for J/ijj can be taken to 
be equal to that of the p meson, as it was done in p (vector meson distribution amplitudes 
were elaborated in [22?]), although this assumption is not completely justified. However, 
we cannot say much about the B meson distribution amplitude, except that it can be 
modeled or extracted form the experimental data j^J, which is again model dependent. 
Fortunately, after some simplification, the result depend only on the first moment of the 
4>b, Jq d^(pB(0 1 '£ = m B I ^Bi and therefore there is a need for fixing just one parameter, A#- 
There is not much known about this parameter, except its upper bound, 3A# < 4(mB — rrib), 
or effectively, X B < 600 MeV \l4\. 

Here, we would like to discuss our results for the soft nonfactorizable contributions in 
comparison with the hard nonfactorizable effects calculated in QCD factorization approach. 
As it was already noted in jlQ| . in the heavy quark limit the soft nonfactorizable contribu- 
tions are suppressed by 1/mb in comparison to the twist-2 factorizable part, which confirms 
the expansion in ()70|). With the inclusion of the hard nonfactorizable corrections, the ci2 
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parameter (JT5j) appears as follows 



F+(m 2 



a s F hard {^) + 



(71) 



The hard nonfactorizable contribution F hard were calculated in [5J]. The analysis was 
done up to twist-3 terms for the K meson wave function which enters the calculation of the 
T 11 hard scattering kernel in (fTT)|) . It is a well known feature of QCD factorization that it 
breaks down by inclusion of higher-twist effects. The hard scattering kernel T 11 becomes 
logarithmically divergent, which signalizes that it is dominated by the soft gluon exchange 
between the constituents of the J/ip and the spectator quark in the B meson. In the QCD 
factorization this logarithmic divergence is usually parameterized by some arbitrary complex 
parameter r as J^dv/v = \n{m B / Aq C d) + r and although it is suppressed by l/m b , this 
contribution is chirally enhanced by a factor 2m 2 K / ((m s + m u )). This large correction makes 
it dangerous to take the estimation for the twist-3 contribution literally, due to the possible 
large uncertainties which the parameter r bears with. 

The estimation done in the QCD factorization [5J shows hard-gluon exchange corrections 
to the naive factorization result of the order of ~ 25%, predicted by the LO calculation with 
the twist-2 kaon distribution amplitude. Unlikely large corrections are obtained by inclusion 
of the twist-3 kaon distribution amplitude. Anyhow, due to the obvious dominance of soft 
contributions to the twist-3 part of the hard corrections in the QCD factorization |l(, it is 
very likely that some double counting of soft effects could appear if we naively compare the 
results. Therefore, taking only the twist-2 hard nonfactorizable corrections from [5] into 
account, recalculated at the \ih scale, our prediction changes to 

|« 2 | = 0.17 — 0.19 | M=Aib (72) 

The prediction still remains to be too small to explain the data. 

Nevertheless, there are several things which have to be stressed here. Soft nonfactorizable 
contributions are at least equally important as nonfactorizable contributions from the hard- 
gluon exchange, if not even the dominant ones. Soft nonfactorizable contributions are of 
the positive sign, and the same seems to be valid also for the hard corrections. While hard 
nonfactorizable corrections have an imaginary part, in the calculation of soft contributions 
the penguin topologies as potential sources for the appearance of an imaginary phase were 
not discussed, but they are expected to be small. 
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VIII. CONCLUSION 



We have discussed the nonfactorizable contributions to the B — > J/ipK decay and have 
calculated leading soft nonfactorizable corrections using QCD light-cone sum rules. In spite 
of theoretical uncertainties involved by application of LCSR method to the B — > J/ipK 
decay discussed in Section El and a possible influence higher charmonium resonances to the 
sum rule, the predicted correction clearly favors the positive value for Fg K and therefore of 
a 2 . 

Recent first observations of the color-suppressed decays of the type 5° — > D^°ir° by 
CLEO Q and BELLE Q indicate also the positive value for a 2 parameter. Although 
these data show that a 2 is a process dependent quantity, which is clearly exhibited by the 
difference in the prediction for a 2 in B° — > D^ ir° and B — > J/ipK decays by almost a 
factor 2 (\a 2 (B° -> D^h°)\ = 0.57 ±0.06 vs \a 2 (B -> J/^)| = 0.28 ±0.03), the positive 
value for a 2 can be clearly deduced in both cases. This is just opposite to the predicted 
negative values of this parameter in D meson decays. The tendency to a positive value of 
a 2 in B decays was also observed in the global fit of decay amplitudes to the data where 
the arguments in favor of a sign change of a 2 from negative to the positive when going from 
D to B decays were presented. 

Moreover, these recent experimental results on 5° — ► D^n point out large nonfactoriz- 
able contributions, as well as the large final state interaction phases in the color-suppressed 
(class-II) decays j22|. Soft corrections obtained in this paper add up to this picture, being 
significantly larger than soft corrections in the B — > tttt decay. 
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